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Abstract 

A new formulation is derived for the commutator-errors in large-eddy simula- 
tion of incompressible flow. These commutator-errors arise from the applica- 
tion of non-uniform filters to the Navier-Stokes equations. As a consequence, 
the filtered velocity field is no longer solenoidal. The order of magnitude of the 
commutator-errors is compared with the divergence of the turbulent stresses. 
This shows that one can not reduce the size of the commutator-errors inde- 
pendently of the turbulent stress terms by some judicious construction of the 
filter operator. Similarity modeling for the commutator-errors is presented, 
including an extension of Bardina's approach and the application of Leray 
regularization. The performance of the commutator-error parameterization is 
illustrated with the one-dimensional Burgers equation. For large filter-width 
variations the Leray approach is shown to capture the filtered flow with better 
accuracy than is possible with Bardina's approach. 
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The desire to extend large-eddy simulation (LES) to complex flows generally implies that 
one is confronted with strongly varying turbulence intensities within the flow-domain and 
also as a function of time. In certain regions of the flow a nearly laminar flow may arise while 
a lively, fine-scale turbulent flow can be present simultaneously in another region. In the 
filtering approach this can be accommodated using a filter operator with non-uniform filter- 
width that may depend on both space and time. The use of such filters, however, further 
complicates the subgrid closure problem through the appearance of additional commutator- 
errors [1] . We will formulate a systematic modeling of the dynamics of these contributions. 

Distinguishing between which flow-features are 'subgrid' and which are 'resolved', de- 
pends on the local filter- width A. Spatial and temporal variations in A therefore imply 
additional energy transfer mechanisms among the scales in the flow, besides the well-known 
energy-transfer due to the quadratic nonlinearity in the Navier-Stokes equations. If a flow 
structure propagates from a region of small filter- width into a region with strongly increased 
filter-width, it would appear as if part of this structure would turn from a 'resolved' to a 
'subgrid' feature, merely by translation. The reverse can also be imagined, leading to the 
apparent emergence of resolved structures from a collection of subgrid scales. This sug- 
gests additional sources of local energy drain or backscatter, depending on the specific local 
filter-width variations in the direction of the instantaneous local flow which require explicit 
parameterization. 

The traditional use of convolution filters in LES necessarily imphes that the width of 
the filter is constant. However, efficient extension of the LES approach to turbulent flows in 
complex geometries or to cases with strong spatial variation of turbulence intensities, calls 
for the introduction of non-uniform filter-widths. The starting point is an extended filter L 
which, in one spatial dimension, is defined by 



where A± denote the upper and lower bounds and A = A_|_+A_. In complex flow geometries 
the variations in turbulence intensities poses different local requirements on the amount of 




(1) 
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detail with which the flow should be represented. Such a situation can be formulated by 
allowing a non- uniform filter- width as given in (1). The application of such filters gives rise 
to a number of extra closure terms in addition to the well-known turbulent stresses. 

If one applies the filter (1) to the incompressible fiow equations, commutator-errors arise 
since dxU ^ d^u or, written differently, L{dxu) — dx{L{u)) = [L,dx\{u) 7^ in terms of the 
commutator of L and the derivative operator Ox- For the filtered continuity equation we 
find 

djUj^-[L,dj\{uj) (2) 

Hence, the filtered continuity equation is no longer in local conservation form and variations 
in the filter-width imply that Uj is not solenoidal. Filtering the Navier-Stokes equations in 
the same way yields 

dtUi dj{uiUj) + diP - -^djjUi = -{[L, dt]{ui) 
+ dji[L,S]{ui,Uj)) + [L,dj]{S{ui,Uj)) 

+ [L,d,]{p)-^[L,d,,]{u,)) (3) 

We observe that commutators emerge of filtering and the product operator S{f,g) = fg 
as well as commutators of filtering and first and second order partial derivatives. Filter- 
ing a linear term such as dtUi gives rise to a 'mean-fiow' term dfUi and a corresponding 
commutator-error [L,dt]{ui). Filtering the nonlinear convective terms leads to the diver- 
gence of the turbulent stress tensor = [L, S]{ui,Uj) and an associated 
commutator-error [L, dj]{S{ui, Uj)). The local conservation form of the Navier-Stokes equa- 
tions is no longer maintained, in the same way as observed in (2). 

The commutators in (2) and (3) satisfy algebraic identities. If we consider any two filters 
Li and L2 then 

[LiL2,5] = [Li,5]L2 + Li[L2,S'] (4) 

which is known as Germano's identity [2]. Likewise, 
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[Li, [L2, S]] + [L2, [S, L^]] + [S, [Li, L^]] = (5) 

which is interpreted as Jacobi's identity. These identities are also satisfied by [L, dt] and 
[L, dj] which shows that the structure of the LES closure problem is closely related to the 
Poisson-bracket in classical mechanics. In that context Germano's identity is known as 
Leibniz' rule. The identities (4) and (5) can be used to guide (dynamic) subgrid modehng 
of the central commutators. 

Filtering the incompressible flow equations gives rise to an 'LES-template' in which the 
'Navier- Stokes' operator on the left hand side of (3) acts on the filtered solution {ui,p}. 
In addition, a number of unclosed terms arises of which only the parameterization of the 
turbulent stresses attracted considerable attention in literature. However, the subgrid mod- 
eling problem associated with non-convolution filters entails various additional commutator- 
errors. These terms require explicit modeling in case the spatial and temporal variations 
of the filter-width are sufficiently large. For steady filter-width distributions, to which we 
restrict ourselves here, the magnitude of these contributions can be quantified in terms of 
u • VA = UjdjA. 

The dynamic effects of the commutator-errors have been considered unimportant by some 
authors, provided a suitable class of filters would be adopted. In [3] such a class of filters was 
constructed and the commutator-errors corresponding to these filters could be made of high 
order in A. Likewise, [4] considers the commutator-errors to be of minor importance in case 
high order filters are used. Although it is correct that the commutator-errors can be made 
small with the proper filter, one has to realize that with the same filter the divergence of the 
turbulent stress tensor is also reduced to the same order in A. The use of higher order filters 
would hence only imply a gradual convergence to the unfiltered Navier-Stokes equations. It 
is not possible to reduce the size of the commutator-errors independently of dj{[L, S]{ui, Uj)) 
merely by constructing suitable filters. The subgrid modeling of the dynamic significance of 
the commutator-errors therefore remains a largely open problem. 

In order to establish the importance of the commutator-errors relative to the turbulent 
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stress contributions we introduce general A'"-th order filters by requiring L{x^) — x"^ for 
m — Q, 1, . . . , N — 1 [4]. Application of such a filter yields: 

U{X)-U{X)^ {A"'{x)Mrr,{x))u^^\x) 
m>N 

where u^"^^ denotes the m-th derivative and Mm{x) is related to the m-th moment of L. To 
leading order u — u ^ . For the commutator-error we find 

[L,d,]{u) = - J2 (A™MJ'«(-) 

m>N 

where the prime denotes differentiation with respect to x. The commutator [L, S]{u) can be 
expressed as: 

[L, S]{u) = ^ (A-M„) _ 2„^M) _ _ 

m>N 

The scaling of the turbulent stresses with is readily verified for > 1. In case N = 1 
the commutator scales with A^ since (u^)' — 2uu' . Combining these expressions one may 
obtain the leading order behavior of the flux terms for symmetric filters as: 

[L,d^]{S{u)) ^ A{x) (A^-^A'M^^) (A^M^) 
d,{[L, S]{u)) ^ a{x) (A^-^A'Mjv) + h{x) (A'^M'j,) 

+ c{x) (A^M^) (6) 

where A, B, a, b and c are smooth, bounded functions which contain combinations of 
derivatives of the solution u. For non-uniform filters [L,dx]{S) ~ A'A"^"^ which is the 
same leading order behavior as for dx[L,S]. Hence, it is not possible to remove only the 
commutators [L, d^] by a careful construction of the filter [3]. In fact, all filters that would 
reduce [L, dx] are of higher order and consequently will also reduce [L, S] with the same 
order. 

One may use a Fourier-mode analysis to relate the dynamic significance of the 
commutator-errors to variations in A and the wavenumber k of the mode. For second 
order filters such as the top-hat or Gaussian filter one has 
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\\[L,d.]{S{um \A'/A\ 

\mL,s]{um^ \k\ 

where || • || denotes the L2-norm. This shows that if |A'| <S |A;A| then filter- width non- 
uniformity can be disregarded and it should be sufficient to model only Tij. This shows that 
only strongly bounded variations in the filter-width will reduce the size of the commutator- 
errors significantly while keeping the magnitude of the turbulent stresses unaffected. If, 
however, for efficiency reasons or due to, e.g., wall-proximity, a sufficiently smooth variation 
of A is not possible, one has to resort to direct modeling of the commutator-errors. 

In the absence of a comprehensive theory of turbulence and its non-uniform spatial 
and temporal representations, the modeling of the turbulent stresses and the commutator- 
errors relies to some degree on limited empirical knowledge. Here we restrict ourselves 
to similarity modeling and consider two different approaches. Specifically, we will extend 
Bardina's approach [5] to include commutator-errors and wc derive the implied subgrid 
models arising from Leray regularization [6]. 

Bardina's similarity model for the turbulent stress tensor arises from applying the defi- 
nition of Tij to Ui, i.e., 

Tij mfj = [L, S] (ui, Uj) = UiUj - (8) 

Extending this idea to the commutator-error suggests the following parameterization: 
[L,dj\{uiUj) [L,dj\{uiUj). In [7] this model showed a high correlation for turbulent 
boundary layer flow which partially substantiates this approach. 

Recently, the Leray regularization principle [6] was revived in the context of LES [8] . In 
this approach the convective fluxes are replaced by UjdjUi, i.e., the solution u is convected 
with a smoothed velocity u. The governing Leray equations can be written as [6] 

djUj = ; dtUi + UjdjUi + diP — -^djjUi = (9) 

Re 

This formulation can be written in terms of {ui,p} in case we assume a (formal) inverse 

of L, i.e., Uj = L^^iuj). After some calculation one obtains the filtered continuity equation 

(2) and the filtered momentum equation as 
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dtUi + djiuiUj) + dip - —djjUi = -{[L, dt]{ui) 



+ {djirriij) + UidjUj] + [L, dj]{S{ui, Uj)) 
+ [L,a,](p)-i-[L,a,,](«,)) 



(10) 



The divergence of the turbulent stress tensor in (3) is represented in terms of the asym- 
metric, filtered similarity-type Leray model m^- and an additional term associated with the 
divergence of the filtered velocity field: 



where mfj — UjUi — UjUi [8] and the commutator-error is expressed as [L,dj\{uiUj) — > 
[L^dj\{jljUi}. The other commutator-errors are identical to those in (3) with the under- 
standing that in actual simulations every occurrence of an unfiltered flow-variable implies 
the apphcation of to the available field. The Leray model was shown to provide good 
predictions of three-dimensional turbulent mixing at arbitrarily high Reynolds number using 
a uniform filter [8]. 

To assess the effects of the commutator-errors and determine the quality of the Bardina 
and Leray modeling we consider the one-dimensional Burgers equation. This provides a 
model-system which has the same basic structure under filtering as expressed in (3). All 
relevant commutators appear in the filtered Burgers equation. The initial solution is a 
Gaussian profile which rapidly develops into the well-known 'ramp-cliff' structure. We use 
Re = 500 to obtain a sharply localized cliff region, and apply periodic boundary conditions. 
Explicit time-integration, restricted by stability time-steps, and second order accurate spatial 
discretization are adopted. To avoid numerical errors we use high spatial resolution, typically 
with N = 2048 intervals. Explicit filtering is done with trapezoidal quadrature applied to 
the top-hat filter. 

We consider a non-uniform grid with grid-spacing hi — {i/N)(l+gi) where £ is the length 
of the domain. The grid is chosen to be non-uniform only in an interval around i = N/2 in 
computational space. For the illustrations we use 



(11) 
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Qi — A sin ( 271 



(i - N/2) 



) 




(m - q)N 



{N{m - 2q)/m) 



m 



and otherwise. Since I] fl'i = this grid preserves the end-points. The parameters q and m 
control the region where the grid is non-uniform. Here, we use g = 3 and m = 8. The local 
filter- width Aj = — where n is chosen such that /S. — 1/D with = 8 or 16 in the 
uniform regions of the grid. With N = 2048 this implies n = 128 or 64 respectively. The 
parameter A <1 controls the ratio between largest and smallest intervals {1 + A)/{1 — A). 

In figure 1 we collected the contributions to the total convective flux for a representative 
uniform and non-uniform case. We decomposed the convective flux as 



identifying on the right hand side the 'mean' flux, the 'SGS-fiux' and the 'commutator-flux' 
respectively. In figure 1 the solution and the filtered solution are included displaying the 
'ramp-cliff' structure. The total flux in figure 1(a) is piecewise linear and the SGS flux is 
localized in the cliff-region where filtering is effective. In figure 1(b) there is a significant 
distortion of the filtered solution due to the filter-width non-uniformity. Two characteristic 
contributions due to the commutator-error arise. On the 'ramp side', the non- uniform filter- 
width near x — —?> strongly influences the mean flux. The commutator-error compensates 
for this such that the total flux remains nearly linear in x. Within the 'cliff-region' the 
commutator-flux is comparable to the SGS-flux. 

In figure 2(a) we show the locations of the front and back of the ramp-cliff solution as 
a function of time. These locations are defined where \u\ equals £max(|M|) with e = 0.05. 
Comparing filtered Burgers results with predictions from the Leray and Bardina parame- 
terizations, the Leray results are more accurate. This was confirmed by considering the 
minimal and maximal values of u which are also better predicted by the Leray model. The 
L2-norm of the fluxes for these cases show that the commutator-flux is about 1/3-1/2 the 
value of the SGS-flux. In cases with smaller grid non-uniformities the direct Leray modeling 
still enhances the accuracy of the predictions notably. The Leray model is considerably less 
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expensive than the Bardina model and it also better preserves qualitative properties of the 
filtered Burgers solution cf. figure 2(b). The Bardina parameterization creates additional 
structure in the solution, which is not present in the filtered Burgers result. 

The commutator-errors have been expressed as commutators of filtering and partial 
derivatives. The magnitude of the commutator-errors can not be reduced independently 
of the SGS-fluxes. Instead, for sufficiently large grid non-uniformities direct modeling is 
needed. For the one-dimensional Burgers equation the Leray parameterization combines 
computational efficiency with high accuracy. This motivates the use of the Leray commu- 
tator model in more complex flows. The a priori specification of the spatial filter-width 
variations is not generally possible for complex cases. Therefore, the local filter-width needs 
to be related to the resolved solution to facilitate a dynamic response of local filtering and 
the evolving flow. This is subject on ongoing research. 
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FIGURES 

FIG. 1. Snapshot of the solution (multiphed by 1/2) (soUd) and filtered solution (solid; markers 
o). Convective flux: total (dots), mean (dash-dotted), turbulent stress (dashed), commutator-error 
(solid with *). In (a) we use A = i/lQ and in (b) the non-uniform case with A = 1/2 is shown. 
Underneath in (b), the grid-spacing (minus 0.2) as a function of x is presented. 

FIG. 2. Location of the head of the cliff (upper curves) and the tail of the ramp (lower curves) 

in (a) and in (b) snapshot of the filtered solution: filtered Burgers (solid), Leray (dashed) and 
Bardina (dash-dotted) for A = 0.85 
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FIG. 1. Geurts, Holm: Phys.Fluids 
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